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Riometer observations of ionospheric absorption at 10 and 15 Mc/s in both the ordinary and 
extraordinary modes were conducted in Alaska during the eclipse of July 20, 1963. The equipment 
was located so that the path of totality at 70 km height was directly overhead, and considerable reduc- 
tion was detected from the ambient daytime absorption, which fell to some 40 percent of its normal 
value. 

The behavior of the D region during an eclipse is considered theoretically, and the observations 
are analyzed in terms of this treatment. It is shown that the observations can be explained in terms 
of ionospheric parameters which, with one exception, accord with values given by other methods. 
The use of future eclipses in Z)-region studies is briefly discussed. 



1. Introduction 

During the solar eclipse of July 20, 1963, riometers 
on frequencies of 10 Mc/s and 15 Mc/s were operated 
near the path of totality at the Fishhook Bar, near 
Palmer, Alaska (61° 41.5' N, 149° 14.3' W, geogr.). 
The riometers were dual-polarized, recording both the 
ordinary and extraordinary components of the cosmic 
radio noise; the antenna beams were symmetrical 
about a vertical axis with maximum response at the 
zenith and half-power gain at 30° off axis. Using the 
methods of Little, Lerfald, and Parthasarathy [1964], 
it is possible from such records to determine the 
absorption suffered by each component in traversing 
the terrestrial ionosphere. 

The absorption in decibels of a radio signal of effec- 
tive angular frequency av, traversing a path through 
the ionosphere, is given by 



A(a> e )- 



Jo 



A(cue, h) N{h) dh 



(1) 



where k((*) e , h) is the absorption per unit path length 
effective at angular frequency w e for unit electron 
density, and N(h) is the electron density as a function 
of height. The absorption efficiency k(a) e , h) is a 
complex function of the electron collision frequency 
v(h), the strength of the terrestrial magnetic field, 
and the angle between the magnetic field vector and 
the direction of propagation of the radio wave. The 
form of the magneto-ionic equations given by Sen and 
Wyller [1960] have been used to compute absorption 
values in this paper. Figure 1 shows the absorption 
efficiency values as a function of height for the case of 
radio waves penetrating the ionosphere vertically. 
The effective angular frequency a> e = co±a>L, where a> 
is the receiving frequency, o>l is the longitudinal com- 
ponent of the angular gyrofrequency and the + and — 



pertain respectively to the ordinary and extraordinary 
polarizations. Because of the variation of absorption 
efficiency and electron density with height, the absorp- 
tion occurs mainly in the height range 60 to 90 km. 
For this reason the riometers were situated somewhat 
south of ground-level totality so that the eclipse would 
be total in the D region overhead. 

Thus, the riometer records a quantity closely related 
to the average electron density in the D region. During 
the eclipse the intensity of ionizing radiation is re- 
duced, and the normal processes of recombination 
bring about a reduction of electron density. The 
object of the experiment was to investigate what can 
be learned from riometer data about these recombina- 
tion processes. Previous investigations of recombina- 
tion in the D region have, in any case, been mainly 
theoretical; experimental evidence on the subject is 
badly needed. 

Absorption measurements by riometer had pre- 
viously been made by Gordon [1961] during the eclipse 
of October 2, 1959. Observing conditions were 
unfavorable at the time, since totality occurred just 




| -IO iq-13 , -M , -l„ | -.* 

ABSORPTION EFFICIENCY IN dB/m /(ELECTRON/m 3 ) 

Figure 1. Absorption efficiency k(o>, h) for ordinary and extraor- 
dinary signals at 5, 10 and 20 Mc/s, and for the mixed mode at 
30 and 50 Mc/s. 
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at sunrise; he used a 30 Mc/s riometer with an antenna 
directed toward the horizon. This antenna orientation 
resulted in several times more absorption than could 
be expected on a vertically directed antenna, but made 
quantitative interpretation very difficult. It is obvious 
that results from equipments using zenithal antennas 
of known polar diagram are more amenable to inter- 
pretation. The use of two frequencies, moreover, 
opened the possibility of finding the distribution of 
electron density with height in this region [Partha- 
sarathy, Lerfald, and Little, 1963]. Conditions were 
still not ideal, however, because the eclipse path was 
near the auroral zone, with a consequent risk of con- 
tamination by sporadic ionization from auroral par- 
ticles. 

We shall first discuss the behavior to be expected 
during an eclipse, and will then present the observa- 
tional results from the eclipse of July 20 and discuss 
their interpretation. 

2. The Prospect: Theory of D-Region Be- 
havior During an Eclipse 

The following processes will be taken into account. 

(i) Production of electrons by ionizing radiation: 

The total rate of production is the sum of contri- 
butions from electromagnetic and corpuscular emis- 
sions: the former from the visible solar disk and per- 
haps also from the corona (source of soft x rays); 
the latter from the steady flux of cosmic rays and, when 
the observations are made at high latitude, from 
sporadic auroral particles. 

(ii) Formation of negative ions: 

Negative ions are formed by the attachment of elec- 
trons to neutral particles, and are dissociated by the 
processes of photodetachment and collisional detach- 
ment. In the equilibrium state the density ratio of 
negative ions to electrons is 



k = 



p + Kn 



(2) 



where a is the attachment coefficient, p the photo- 
detachment coefficient (which depends on the sun's 
intensity in the visible and near infrared), k the 
collisional detachment coefficient, and n the number 
density of neutral molecules. 

(hi) Recombination with positive ions: 

Both the electrons and the negative ions can be 

removed by recombination with positive ions. If the 

rates of these processes are respectively a e and w, 

the "effective recombination rate" of the electrons is 



a = a e -\- ka^. 



(3) 



In terms of these processes, the changes of electron 
density can readily be shown [Nicolet and Aikin, 
1960] to be described by 



d/V. 

dt ' (1 + X) 



._ Q 



-aN 2 



N dk 
"(1 + X) dt 



The derivation of (4) assumes that the electrons and 
negative ions come to equilibrium rapidly compared 
with changes of q: that is, that (2) applies throughout. 
As the time constant of process (ii) is only about one 
second in the D region, the assumption is valid for the 
present analysis. 

2.1. Primary Effects of the Eclipse 

The first effect of the eclipse is to modify the rate 
of production of ionization, q. Only the electromag- 
netic component is changed, so we may put 



q = Fq + q c 



(5) 



where q is the electromagnetic component from the 
unobscured sun, F is the fraction of the sun exposed, 
and q c is the corpuscular component which, not coming 
from the sun, is unaffected by the eclipse. If the 
solar component is confined to the visible disk and 
uniformly distributed across it, F changes between 
unity and zero during the eclipse, as shown in the case 

— = of figure 2. (The unit of "eclipse time" in fig. 2 

9° 

is the time for the moon to move by one radius relative 

to the sun.) If any appreciable contribution comes 
from the corona, then the effective solar radius will 
be slightly increased and F will be slightly different. 
Also, active regions on the sun could result in irregu- 
larities in the change of solar ionizing flux during the 
course of the eclipse. At low latitude the corpuscular 
component q c is probably constant during the eclipse, 
but this cannot be assumed at high latitude unless 
there is evidence that auroral activity is absent. 

The second effect of the eclipse is to modify the 
photodetachment coefficient, and hence the ratio of 
negative ion density to electron density, k. Thus 



k = 



Fpo + K> 



(6) 



where p is the photodetachment coefficient for unob- 
scured sun and F the exposed fraction of the solar 
disk. For convenience we have put an 2 =A and 
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FIGURE 2. Production function, Q. 

( 4) For the eclipse of July 20, 1963, one eclipse unit == 36V2 min., duration of eclipse = 146 min, 

and duration of totality = 93 sec. If q c — 0, Q is just the exposed fraction of the solar disk, F. 
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kti = K. During the eclipse, therefore, 
between day and night values: 

A - A X -±. 



At a sufficiently great height, X is always negligibly 
small. In the D region, however, it is probably signifi- 
cant at all times. For a height of 70 km, for instance, 
Nicolet and Aikin [1960] give X D = 0.6; Bailey [1959] 
estimates Xd — 40, X^ = 300. 

It seems therefore, that both the recombination and 
the attachment processes may well affect the obser- 
vations. Nevertheless, we will first treat them sepa- 
rately so as to study the peculiar effects of each 
process. 



2.2. Behavior if X is Constant During the Eclipse 
If k = kp throughout the eclipse, (4) becomes 
dN_ q 



k changes days; and it may be possible to estimate the ratio 



qc 



dt (l + \z>) 



■a D N\ 



(7) 



kp and an denoting the daytime values. The equi- 
librium electron density by day (No) is given by 



*dN% = 



Qd 



(l + M 



(8) 



where q D is the daytime rate of production. Whence, 
using (5) we obtain the convenient form 



dt \N D , 



--a D N l 



°H£J} 



(9) 



where 



F + 



qc 



1+* 



Qo 



Solutions of this equation are given in figure 3 for 

Qc 
different values of the parameter oldNd and for -j- = 0, 

0.1, 0.3, and 1.0. According to Nicolet and Aiken 

[1960], a D at 70 km is 6 • 10~ 7 cm 3 sec" 1 and N D is 

about 300/cm 3 . Thus, apNp is about 2 • 10~ 4 sec -1 , 

which puts it in the range which should give observable 

effects. If oldNd happens to be too small, say less 

than about 10~ 5 sec -1 , no effect will be seen; and if it 

is too large, say greater than about 10~ 2 sec -1 , the curve 

N 
— will be indistinguishable from the square root of 



the production function. In the intermediate range 
an actual measurement of ap becomes possible if Nd 

and J 2 are known. Nd can be determined from the 
daytime level just before the eclipse or on adjacent 



from the variation of absorption with zenith angle. 

Uncertainty about the value of ~ or about the form 

of F will obviously make the determination of otpNp 
less certain, but the tail of the recovery curve might 
be helpful in this case. After last contact, when 



F + 



qc 
Qo 



= 1, 



1+f 

(7) has the solution 

tanh" 1 (—J - tanh" 1 h^j = a D N D t (10) 

where N = N\ at time t = 0. In this way a determina- 
tion independent of the form of the production func- 
tion can be made — provided, of course, that a recovery 
tail is observed. If a tail is not observed it may only 
be possible to say that otpNp must be greater than or 
less than certain values. 

As riometer observations over a band of frequencies 
can potentially be interpreted as profiles of electron 
density as a function of height [Parthasarathy, Ler- 
fald, and Little, 1963], the variation of an with height 
might also be obtained. The electronic and ionic 
recombination coefficients, a e and a*, could then be 
obtained by means of (2) and (3), making use of the 
strong variation of k with height. 

2.3. Behavior if X Changes During the Eclipse but 
Recombination Is Slow 

Now let us suppose that ap ~ 0, so that there is 
virtually no recombination with positive ions for the 
duration of the eclipse. As the electron density (over 
periods of a second or more) is always related to the 
positive-ion density (/V+) by 



W = 



(1 + X) 



(11) 



and since A^+ is constant in the absence of recombi- 
nation, it follows that 



N = l + k D 
N D 1 + X 



= (1 + Ad) 



1 + 



i+f {£H 



(12) 



Then measurements of — as a function of F will pro- 

N D ^ 

vide values for both kp and X#, from which the ratios 

A K 

— and — may be determined. If the changes of kp and 

X# with height are also known, we expect to find from 

A K 

them values of 5- and -j which vary respectively as the 
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second power and as the first power of the molecular 
density. This would provide some verification of the 
accuracy of the method. 

2.4. Behavior if X Changes During the Eclipse and 
Recombination Is Also Important 

In general, both the recombination and the attach- 
ment processes will occur in the D layer. If the 
changes due to recombination happen to be small, 
then the effects of the two processes may be com- 
bined simply by multiplying each point of the curves 

in figure 3 by the factor D • The modification is 

(1 + X) 

symmetrical about totality. This may also be satis- 
factory as a first order treatment of other cases, but 
the general situation is complicated by the influence 
of X on the effective production and recombination 
rates. This can been seen by writing (4) in terms of 
positive ion density: 



N+ = (l + k)N 



(13) 



The basic recombination process is obviously in- 
fluenced by changes of X unless (it can be shown), 
(a € — otj) is sufficiently small compared with (a e +k(Xi) 
(1H-X), or unless X changes by only a small amount. 
The various coefficients are not well enough estab- 
lished to say whether or not these conditions will be 
met. 

Equations (4) or (13) can be solved numerically, and 
examples of such solutions are shown in figure 4. It 
is seen that a marked dip at totality is introduced into 
a curve which otherwise would pass through a mini- 
mum somewhat later than totality. 

2.5. Some Interpretive Difficulties 

Given an experimental curve showing the variation 

N 
of — during an eclipse, the problem is to isolate three 

kinds of variations: 

(i) The change of ionizing flux during the eclipse; 

(ii) the change of X during the eclipse, which will 

provide information about the attachment and 

detachment coefficients; 

(in) the delayed change, giving information about 

the recombination coefficients. 

If it could be assumed that the only ionizing radia- 
tion was from the visible solar disk and was uniformly 
distributed across it, the analysis would be relatively 
simple. The early and late parts of the eclipse, where 
X is near to its daytime value, would be fitted to a theo- 
retical recombination curve from figure 3a. De- 
partures from this same theoretical curve near totality 
would then be interpreted in terms of changes of X: 
the ratio between the two curves could be taken as 
(1+Az>) 



(1 + X) 



in the first instance, and then the values re- 



fined by comparison with numerical computations. 
The presence of additional sources of ionization which 
remained steady for the duration of the eclipse would 
alter the basic recombination curve but not alter the 
effects of attachment. It would be necessary to esti- 
mate in some way a value for -jj- and then to use the 

appropriate set from figure 3. 

In the presence of sporadic auroral ionization or 
gross nonuniformities in the ionizing flux from different 
regions of the solar disk, no precise solution is possible. 
An uncertainty is introduced whose seriousness can 
however be estimated if the general level of sporadic 
activity or nonuniformity is known. The possibility 
of auroral absorption makes observation at high lati- 
tude somewhat undesirable, but if there is no other 
choice the results could be improved by observing at 
several sites along the eclipse path. Then the 
sporadic effects would appear against different phases 
of the eclipse at the different sites, and only the sys- 
tematic component would be taken as the true effect 
of the eclipse. 

Ideally, multifrequency observations during the 
eclipse will give electron densities as a function of 
both time and height, and the analysis would aim to' 1 
interpret the variation of electron density with time 
at each of several heights. Should the observations not 
be precise enough to give height information, the re- 
sults can provide only "effective parameters' 1 for some 
"effective ionospheric layer." In the present state 
of knowledge of the D region, this could still be 
valuable. However, the possibility that some of the 
absorption is in the E and F regions has to be kept in 
mind. Independent evidence on this point would be 
necessary. 

3. The Reality: Analysis of the Eclipse of 
July 20, 1963 

3.1. Observations 

The absorption observed with the four riometers 
(recording the O- and E-modes at both 10 and 15 Mc/s) 
is shown in figure 5. These observations have been 
corrected for the finite beamwidth of the antenna and 
for mutual leakage between the two circularly polarized 
components. A correction has also been applied to 
the data for the apparent absorption introduced by 
F-region reflection of cosmic noise arriving from angles 
near the horizon (the so-called "aperture" or "iris" 
effect). The absorption normally observed during the 
same period of the day (adjusted to the same solar 
zenith angles) derived from about a month of control 
data recorded before and after the eclipse, is also 
shown. It can be seen at once that the eclipse pro- 
duced a marked reduction of absorption. The records 
also display a fair amount of irregularity, which could 
be caused either by auroral activity which would in- 
crease the absorption, or by interfering signals which 
would apparently decrease it. The former is believed 
to have been the more serious during these observa- 
tions. Whatever its cause, the irregularity presents 
a fundamental limit to the accuracy of any inter- 
pretation. 
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Figure 3. Relative electron density (rr~) at constant X. 

For generality, the parameter adVo is expressed in units of eclipse time. The values are repeated in brackets in sec _1 for the eclipse of July 20, 1963. 
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Figure 4. Relative electron density [t^~j with changing k. 

a D /V,>= 1.37; a e /a, = 5.0. 
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FIGURE 5. Absorption during eclipse of July 20, 1963. 



3.2. Height of Ionization Being Observed 

Because of the irregularity of the records it is not 
possible to work out any changes in the Z)-region pro- 
file during the course of the eclipse. Our interpreta- 
tion will therefore be considered, in the first instance, 
to describe the behavior during the eclipse of an 
absorbing region lying between 60 and 80 km, these 
heights having been derived from the theoretical 
analysis of Lerfald and Parthasarathy (to be published). 
In this analysis the variation of absorption as a func- 
tion of solar zenith angle during periods free from 
auroral activity has been used to investigate the valid- 
ity of various ionospheric models. The N(h) profiles 
between 50 and 90 km resulting from assumed ionizing 
fluxes and rate coefficients were computed using the 
method given by Nicolet and Aikin [1960] and the 
choice of model parameters narrowed by rejecting 
those giving results which are inconsistent with the 
observed variation of absorption with solar zenith 
angle. (The N(h) profiles above 100 km were obtained 
from ionosonde data by the method of Paul and 
Wright [1963].) 

Such an analysis leads to the conclusion that the 
appropriate profile lies between the two given in figure 
6a. The corresponding profiles of absorption against 
height are shown in figure 6b. The absorption is a 
maximum between 70 and 75 km, and exceeds half the 
maximum value throughout the approximate range 60 
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FIGURE 6. Estimated profiles. 

The profile at the time of the eclipse is thought to lie between (i) and (ii). 



to 80 km. There is a smaller E-region contribution 
above 100 km, which may amount to a tenth of the 
Z)-region contribution. In view of other uncertainties, 
no attempt was made to correct for this ^-region 
absorption contribution, although it would be possible 
to do so. 

3.3 Interpretation 

From the observations of figure 5 a pattern has been 
abstracted to represent the effect of the eclipse. This 
eclipse pattern was derived by dividing the absorption 
on the eclipse day by the corresponding absorption on 
normal days, and taking an average over the four ob- 
serving frequencies. The average was biased towards 
the channel showing the smallest absorption, and ef- 
forts were made to smooth out changes probably 
caused by sporadic auroral absorption. The resultant 
pattern is given in figure 7a; this pattern we now wish 
to interpret in terms of the recombination and attach- 
ment processes. 

The recovery of ionization during the latter part of 
the eclipse is relatively slow, and the electron density 
is still well below the daytime value at last contact 
(eclipse time = + 2). This means that recombination 
can be neither very slow nor very rapid. Before fitting 

a curve from figure 3, a value of ^ must be chosen. 

The analysis which led to the estimated profiles of 

figure 6 also indicates that a value of jp between 0.05 

and 0.5 would be most appropriate. The curves for 



: 0.1 and 0.3 have therefore been used, and of these 

the values cxdNd = 0.375 and 0.500 (eclipse units) -1 
give the best fits respectively. These values are 
equivalent to 1.7 • 10~ 4 sec -1 and 2.3- 10~ 4 sec -1 , which 
are somewhat greater than, but still compatible with, 
az) = 6*10~ 7 cm 3 sec" 1 given by Nicolet and Aikin 
[1960] and Nd = 700 cm -3 indicated for 70 km in figure 
6a. It is of some interest that it becomes difficult to 
fit a curve at all if ^ > 0.3; of our two values, ^ = 0.1 

might therefore be preferred. 

As shown in figure 7b, the ratio between the observed 
eclipse pattern and the theoretical one assuming con- 
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stant negative ion-electron ratio (X = X D ) is symmetrical 
about totality. In accordance with the foregoing dis- 
cussion we shall therefore regard this ratio as an indi- 
cation of • Assuming that the solar illumina- 

(1 + Ad) 

tion factor F is as in figure 2, we thus obtain the 

variation of n . with F shown in figure 8. 

(1 -h Kd) 

From the variation of /n , with F may be deter- 

(1 -h Kd) 

mined the day and night values of X: i.e., X# and X#. 
It follows from (12) that 

{ ^ +( «- 1)} [ / ,{i-i-Ui-]= 1 , 



where R = 



1 + A 



To determine Xd and X/v from the 



1 + Xd 

observed variation of R with F, we differentiate this 
equation with respect to F and rearrange it to obtain 



R 



dF) k N - 



dR 



1 



\ D dF 1 + \ D 



(14) 
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FIGURE 7. Eclipse absorption pattern. 

(a) Relative electron density: observed, and theoretical at constant k. 
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R + F — J with 



dR 
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Values obtained from figure 8 are plotted in this form 
in figure 9. The relation is seen to be reasonably 
linear: the intercept gives Xd = 0.35 and the slope 
gives Xat= 1.7. 

4. Discussion 

4.1. Comparison With Other Results 

If we accept that the measurements refer to the 
height range 60 to 80 km, then our estimate of X# agrees 
reasonably well with that of Nicolet and Aikin, who 
give Xz) = 0.35 at 72 km. The estimate of X#, however, 
is smaller than any given previously for this height. 
X^ is a parameter about which relatively little is 
known, mainly because of ignorance of the collisional 
detachment coefficient, k, [see Reid, 1964, for exam- 
ple]. Unfortunately, our derivation is also the least 
accurate part of the present analysis because of the 
possible contribution to the absorption from ionization 
at the higher levels, say above 80 km. Such a con- 
tribution might amount to 20 percent of the total and, 
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being in a region where X is much smaller than at 70 
km, it would tend to persist throughout totality. 
Nevertheless, there is the experimental fact that nearly 
40 percent of the absorption remained at totality; even 
assuming that half this was due to a high-level con- 
tribution, it is hard to make our effective value of A# 
greater than about 10. 

4.2. Value of Riometer Observations During Future 
Eclipses 

It appears that eclipse observations can reasonably 
be applied to the determination of certain ionospheric 
parameters, and such efforts would seem worthwhile 
if only to derive values of X^. This would be particu- 
larly so if changes in the Z)-region electron density pro- 
file against height could be measured. 

The interpretation of absorption measurements is 
admittedly hampered by many side effects such as the 
contributions by the higher ionospheric layers, changes 
in the corpuscular ionizing fluxes, the contribution of 
the solar corona, or nonuniformities of emission across 
the solar disk. However, alternative methods of 
studying the D region used to date have similar 
problems and limitations. Since the change in in- 
tegrated absorption during an eclipse at a frequency 
such as 30 Mc/s is only of the order of 0.1 to 0.2 dB, it 
is not possible to obtain definitive results using con- 
ventional riometers, which typically have a resolution 
of the order of 0.1 dB. Two alternatives are open: 
the development of equipments operating at frequen- 
cies such as 30 Mc/s which can resolve the absorption 
to about 0.01 dB, or absorption measurements at lower 
frequencies to take advantage of the approximate 
inverse-frequency-squared law of absorption. The 
latter choice has the advantage that the total absorption 
change expected and the difference between absorp- 
tion on the ordinary and extraordinary polarizations 
is relatively large, but the method is limited by the 
effects of the higher ionospheric layers and the effects 
of propagated interference. It might be added that 
recent tests of equipment have shown that the con- 
ventional servo type riometer, which uses electronic 
(diode) switches operating at a few hundred cycles per 
second at the input, is more susceptible to interference 
than a nonswitching type of equipment. The latter 
would therefore be preferable for eclipse work. 

The optimum conditions for an experiment of this 
type include low values of E- and F-region critical 
frequencies, an absence of sporadic corpuscular bom- 
bardment and an eclipse with a long duration of totality. 



A geographic location somewhat equator-ward of the 
auroral zone during a sunspot minimum period would 
probably provide the best conditions in practice. 

5. Conclusions 

(1) The eclipse produced notable effects in both the 
10 and 15 Mc/s riometer observations, the absorption 
falling to 0.4 of its usual value for that time of day. 

(2) The effect can be satisfactorily interpreted in 
terms of recombination and attachment processes in 
the D region. The values deduced for the recom- 
bination coefficient and for the concentration ratio 
of negative ions to electrons under day conditions are 
not inconsistent with those determined by other means. 
The ratio determined for night conditions is con- 
siderably smaller than other suggested values. 
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